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1 ABSTRACT 

2 Cell navigation is directed by inhomogeneous distributions of extracellular cues. It is well 

3 known that noise plays a key role in biology and is present in naturally occurring 

4 gradients at the micro- and nanoscale, yet it has not been studied with gradients in vitro. 

5 Here, we introduce novel algorithms to produce ordered and random gradients of 

6 discrete nanodots - called digital nanodot gradients (DNGs) - according to monotonic 

7 and non-monotonic density functions. The algorithms generate continuous DNGs, with 

8 dot spacing changing in two dimensions along the gradient direction according to 

9 arbitrary mathematical functions, with densities ranging from 0.02% to 44.44%. The 

10 random gradient algorithm compensates for random nanodot overlap, and the 

11 randomness and spatial homogeneity of the DNGs were confirmed with Ripley's K 

12 function. An array of 100 DNGs, each 400 * 400 pm 2 , comprising a total of 57 million 

13 200 x 200 nm 2 dots was designed and patterned into silicon using electron-beam 

14 lithography, then patterned as fluorescently labeled IgGs on glass using lift-off 

15 nanocontact printing. DNGs will facilitate the study of the effects of noise and 

16 randomness at the micro- and nanoscales on cell migration and growth. 
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1 Introduction 

2 Gradients are fundamental to many phenomena of biology, from directing axonal 

3 navigation during neural development to the differentiation of stem cells in response to 

4 an injury [1,2]. Gradients may occur as either (i) free diffusion gradients, (ii) substrate- 

5 bound gradients, or (iii) a combination of both [3]. Free diffusion gradients can be 

6 generated by cell-excreted proteins that diffuse in the extracellular space. Substrate- 

7 bound gradients arise when proteins are bound to the extracellular matrix or cell 

8 membranes [4,5,6]. To better understand cellular haptotaxis - the directed movement of 

9 a cell or axonal growth cone along a gradient of a substrate-bound guidance cue - many 

10 techniques to generate surface-bound gradients in vitro have been developed [7]. Two 

11 classes of concentration gradients exist, namely continuous and digital gradients. 

12 Continuous gradients have protein concentration changing in a constant manner, and 

13 are typically produced by adsorbing molecules from a diffusible gradient. Several 

14 parameters of these gradients have been modulated, including the range and slope, 

15 which was linear, exponential or non-monotonic [8]. For example, linear gradients of the 

16 protein BDNF have been shown to modulate neuron polarization and growth [9] and 

17 linear laminin gradients have been shown to orient rat hippocampal axon specification 

18 based on slope [10]. Secondly, digital gradients which have been recently introduced, 

19 are formed by patterning small dots of protein and varying the density of the dots by 

20 changing their size [11], the spacing between them [12], or both [13]. The advantage of 

21 digital gradients is that they are deterministic, and that the local concentration can be 

22 precisely predicted and does not rely on fluorescence measurements which are prone to 

23 error. Digital gradients have been used to study how retinal ganglion cells identify the 

24 stop zone within graded distributions of repulsive EphrinA5 ligands [13]. However, digital 

25 gradients rarely extend over 1-2 orders of magnitude (OM), whereas it is believed that 

26 the dynamic range of gradients in vivo is between 3-4 OM. To overcome these 
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27 limitations, we previously developed digital nanodot gradients (DNGs), where the 

28 spacing between nanodots (200 nm in diameter) was changed in two dimensions to 

29 produce a dynamic range exceeding 3 OM. These designs were implemented using a 

30 low-cost, lift-off nanocontact printing method to pattern substrate-bound gradients of 

31 proteins and peptides. We employed these patterns for adhesion and migration studies 

32 of C2C12 myoblasts on RGD peptide and netrin-1 gradients, respectively [14]. For these 

33 experiments, gradients of 400 * 400 urn 2 were divided into 53 rectangular boxes of fixed 

34 size and density. With this approach, gradients were non-continuous and had stepwise 

35 changes in density. This is most pronounced in low-density regions where the spacing 

36 between nanodots exceeded the dimensions of the box, requiring larger box sizes and 

37 thus creating large "steps" at the bottom of the gradient. This may prove problematic as 

38 cells may fail to sense a discontinuous gradient if they fall into a region of constant 

39 density. 
40 

41 Noise is ubiquitous in biology [15], and modern patterning technologies can be exploited 

42 to introduce defined amounts of noise and randomness into otherwise regular patterns. 

43 The effect of randomness was evaluated in ordered arrays and repetitive patterns with a 

44 constant average density. In one study, disorder was introduced in arrays of 120 nm 

45 dots spaced 300 nm apart by randomly displacing the dot by up to 50% of the spacing to 

46 avoid overlap [16]. While this approach introduced a controlled amount of noise by 

47 restricting the maximum displacement, it was not random since dots were each 

48 contained within the original grid. Nonetheless, cellular adhesion and stem cell 

49 differentiation of osteoblasts were markedly altered as a result of increasing disorder. In 

50 another study, whole proteome analysis of cells grown on the same disordered patterns 

51 resulted in differential expression of certain proteins in the extracellular signal-regulated 

52 kinase (ERK1/2) pathway [17]. Similarly, controlled amounts of topographical noise on 
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53 nanogratings of 500 nm ridges and grooves has shown to effect PC12 neuronal cell 

54 alignment to the gratings, focal adhesion maturation and directionality [1 8]. 
55 

56 Randomness and noise are also highly relevant to directed cell migration. The 

57 stochasticity of chemo- and haptotaxis has been well studied, and is apparent from the 

58 random-walk like traces of migrating cells [19]. It is well understood that biological 

59 gradients, which appear continuous, are in fact quantized since they are comprised of 

60 individual molecules adsorbed to a surface. The distribution of these molecules is not 

61 deterministic, but stochastic at the nanoscale. The engagement of receptors from 

62 migrating cells with these guidance cues has been modeled within a stochastic 

63 framework [19,20]. Random variations in the gradient also occur at the microscale in vivo 

64 from the local accumulation of chemo- and haptotactic molecules that form concentration 

65 puncta [21]. Gradients formed by the expression of receptors from cells embedded in a 

66 tissue a may become non-monotonic as particular cells over- and under- express a 

67 receptor relative to their position in the gradient [22]. Cells navigating through such a 

68 patchwork of microscale deviations must discriminate against local maxima and minima 

69 to sense the overarching gradient slope. It has been suggested that cells alternate 

70 between periods of sensitization and desensitization and are thereby capable of 

71 maintaining an overall response to long-range gradients [23]. 
72 

73 Random pattern generation is well established, however many patterns are not random 

74 in the mathematical sense, and often implement only an approximation of randomness 

75 to varying degrees depending on the needs of the application. Firstly, computers 

76 typically generate 'pseudo-random' numbers that statistically approximate the properties 

77 of true random numbers using an algorithm based on an initial 'random seed' - generally 

78 a random bit. If this random seed is known, the set of pseudo-random numbers can be 
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79 replicated. However, for most applications, pseudo-random numbers are 

80 indistinguishable from a true random number. 
81 

82 So-called random dot patterns are widely used in the reproduction of images and 

83 patterns using displays or printers that impose digitization and pixilation. Pixilation 

84 implies that pixels are positioned on a regular grid. Adjusting grayscale images to binary 

85 (black and white pixels) through halftoning is achieved by adjusting the ratio of black and 

86 white pixels in a local area, but if pixels are turned on and off according an ordered 

87 pattern, rendering artifacts (such as Moire patterns) arise. By locally randomizing the 

88 pixels turned on or off while matching the overall desired grayscale, it is possible to 

89 avoid such rendering artifacts. These 'random-grid' patterns can be generated using an 

90 array of pseudo-random values and a threshold dictated by the local grey value used to 

91 determine which pixels are black or white (Fig. S1) [24]. With this approach, the value of 

92 each pixel is random, while the position of the pixels are fixed to an ordered grid. These 

93 random-grids have found application in the halftoning of images for display or printing 

94 [25] and the generation of random-dot stereograms for depth perception studies [26]. 
95 

96 Random dot patterns also comprise dot patterns where the position of the dots is not 

97 constrained to a grid as described previously, but is instead randomized within a given 

98 space. However, pseudo-random distributions of the dots, when seen from afar, can 

99 appear noisy and inhomogeneous to the eye. Therefore, algorithms involving so-called 

100 quasi random numbers (also called low discrepancy sequences) are used to create 

101 patterns with a more uniform distribution that appears visually smooth [27]. "Quasi- 

102 random" is a broad term that encompasses numbers with a disordered distribution that 

103 lie between a (pseudo-)random distribution and a regular distribution [28]; a quasi- 

104 random distribution may closely resemble an ordered distribution than a random 
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105 distribution upon close examination. Many different algorithms have been developed to 

106 produce quasi-random numbers. As an example, for liquid crystal displays, quasi- 

107 random dot patterns without overlap are generated by mimicking molecular dynamics in 

108 which two "molecules" cannot occupy the same space [29] or by iteratively repositioning 

109 overlapping dots until overlap is minimized [27]. 
110 

111 Both random-grid and quasi-random patterns have been optimized to produce a 

112 pleasing (macroscopic) image to the eye by introducing a controlled amount of 

113 microscale disorder, avoiding both the viewing artifacts of regular arrays and the 

114 inhomogeneity of pseudo-random dot patterns. Cells however are microscopic and 

115 sense their environment at the micro- and nanoscale. At these scales, random grid or 

116 quasi-random distributions can be significantly different from a random one, and may not 

117 adequately reproduce the effect of a random pattern on cells. It is therefore necessary to 

118 establish an algorithm to produce random dot patterns that preserves randomness at the 

119 micro- and nanoscale to study the effect of random vs. regular dot patterns on cell 

120 navigation. 
121 

122 Here, we introduce continuous DNGs eliminating the stepwise density changes of the 

123 previously reported "step" design [14] with (i) ordered and (ii) pseudo-random positioning 

124 of nanodots, as well as (iii) non-monotonic DNGs that can be implemented using either 

125 (i) or (ii). We discuss various strategies to create noisy DNGs and outline the challenges 

126 in forming random DNGs with accurate slopes, and describe a novel algorithm with 

127 pseudo-random positioning of the dots and compensating for random, but predictable, 

128 dot overlap to achieve the desired coverage. The slope of random DNGs was measured 

129 and compared to the programmed density function, and their randomness verified using 

130 Ripley's K-function. We generated an array of 100 ordered and random 400 * 400 urn 2 
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131 large DNGs made of 200 * 200 nm 2 nanodots, including monotonic and non-monotonic 

132 density curves, with a dynamic ranges spanning from 2.14 to 3.86 OM. Non-monotonic 

133 gradients produced here aim to introduce in a quantitative and repeatable manner 

134 "microscale noise" into surface-bound in vitro gradients. The entire array of 100 DNGs 

135 covers a 35 mm 2 area and is comprised of more than 57 million nanodots. The DNG 

136 array was patterned onto a silicon (Si) wafer by electron-beam (e-beam) lithography, and 

137 was transferred onto glass slides by lift-off nanocontact printing of fluorescently labeled 

138 IgGs [14]. The fidelity of the replication process was evaluated by overlaying the DNG 

139 design with fluorescence microscopy images of the printed IgG proteins. 
140 

141 Materials and Methods 

142 

143 Software & Computers 

144 Algorithms were developed in Matlab R2013a (Natick, MA, USA). Scripts of the ordered 

145 and random gradient algorithms are available for download upon request. A spreadsheet 

146 template (Microsoft Excel 2010) of gradient parameters was imported into Matlab for 

147 processing. The output for each gradient was a text file of coordinates formatted as a 

148 Caltech Intermediate Format (CIF) file. CIFs were imported into CleWin Version 4.1 

149 (WieWeb, MESA Research Institute at the University of Twente and Deltamask, 

150 Netherlands). Gradients were exported from CleWin as Bitmap (BMP) image files for 

151 verification of density using ImageJ 1.47 64-bit (National Institutes of Health, USA) and 

152 Matlab. The 100-gradient array was designed in L-edit (Tanner EDA, Monrovia, CA, 

153 USA). All software was run on a 2013 iMac computer with a 3.4 GHz Intel Core i7 

154 processor with 32 GB 1600 MHz DDR3 memory running Windows 7 through a bootcamp 

155 partition. 
156 
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157 Electron-Beam Lithography 

158 A 4" Si wafer was coated with PMMA resist and the 1 00-gradient array was patterned by 

159 e-beam lithography (VB6 UHR EWF, Vistec, Montreal, QC, Canada), followed by 

160 reactive-ion etching (SystemlOO ICP380, Plasmalab, Everett, WA, USA) 100 nm deep 

161 into the Si wafer. 
162 

163 Stamp Fabrication 

164 After cleaning, the Si wafer was coated with perfluorooctyltriethoxysilane (Sigma-Aldrich, 

165 Oakville, ON, Canada) by vapor phase deposition. An accurate polymer copy of the 

166 wafer was obtained after double replication using polydimethylsiloxane (PDMS) and a 

167 UV-sensitive polyurethane as described in [14]. First, a 6 mm layer of 1:10 PDMS (Dow 

168 Corning, Corning, NY, USA) was poured on the wafer inside a Petri dish, followed by 

169 degassing under vacuum in a desiccator for 10 min. Next, the PDMS was cured in an 

170 oven for 24 h at 60 °C (VWR, Montreal, QC, Canada), then peeled from the wafer. To 

171 remove uncured monomers and other extractables, the PDMS replica was submersed in 

172 70% ethanol for 24 h then baked at 60°C for 4 h. Next, a drop of UV sensitive 

173 polyurethane (Norland Optical Adhesive 63 (NOA 63), Norland Products, Cranbury, NJ) 

174 was applied to the PDMS replica and cured by 600 W of UV light (Uvitron International, 

175 Inc., West Springfield, MA) for 30 s. The PDMS was then removed yielding an NOA 

176 replica of the original Si wafer pattern with 200 nm holes [14]. 
177 

178 Nanocontact Printing 

179 A flat PDMS stamp was cured against a perfluorooctyltriethoxysilane treated flat Si 

180 wafer. Following removal of the extractables as mentioned above, the flat PDMS stamp 

181 was inked with a 10 uL drop of phosphate buffered saline solution (PBS) containing 25 

182 ug/mL of chicken immunoglobulin G (IgG) conjugated to Alexa Fluor 488 (Invitrogen, 
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183 Burlington, ON, Canada). A plasma activated hydrophilic coverslip was then placed on 

184 the drop to spread the solution evenly across the surface of the hydrophobic PDMS 

185 stamp during a 5 min incubation period. After rinsing with PBS and double distilled water 

186 for 15 s each, the inked stamps were briefly dried under a stream of N 2 and immediately 

187 brought into contact with a plasma activated (PlasmaEtch PE-50, PlasmaEtch, Carson 

188 City, NV, USA) NOA master for 5 s. The PDMS and NOA were separated and proteins in 

189 the contact regions were transferred to the NOA. The remaining proteins on the PDMS 

190 were transferred to the final substrate by printing the PDMS stamp for 5 s onto a plasma 

191 activated glass coverslip. 
192 

193 Imaging and Analysis 

194 Images of the original Si master were collected using scanning electron microscopy 

195 (SEM, JEOL, Japan). DNGs of fluorescent IgGs were imaged by fluorescence 

196 microscopy (TE2000 microscope, Nikon, Canada and CoolSNAP HQ 2 camera, 

197 Photometries, USA). Images of the Si wafer were captured with a Panasonic Lumix GH3 

198 DSLR equipped with an Olympus M. Zuiko Digital ED 60 mm macro lens. Dark field 

199 images were captured with an inspection microscope (LV150A microscope and Digital 

200 Sight DS-FM camera, Nikon, Canada). 
201 

202 Results and Discussion 

203 

204 Ordered Gradient Algorithm 

205 Ordered DNGs with continuously changing density were programmed by forming 

206 columns of nanodots with equal vertical spacing while varying the spacing between 

207 columns horizontally (Fig. 1 A). The density of the gradient at any point along the length / 

208 is dictated by an input density function D, and is realized by placing a single nanodot into 
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209 a virtual box to form a unit cell. The dimensions of,- of the unit cell at the /* column of 

210 nanodots are given by the square root of the nanodot area A do t divided by the density 

211 value at the given position (Eq. 1). In this algorithm, the size of the nanodot remains 

212 constant while the dimensions of the unit cell vary for each column of nanodots along the 

213 length. 

214 d , = 4 P**- (1) 

215 The unit cell is largest at low densities and decreases at higher densities, matching the 

216 dimensions of the nanodot at a density equal to one. The unit cell dimensions for each 

217 column are calculated using an iterative approach, starting from the first column of 

218 nanodots at length zero. The position along the length for the next column of nanodots is 

219 given by the cumulative sum of unit cell dimensions (Fig. 1B). 
220 




221 

222 Figure 1: Ordered gradient algorithm for continuously changing density in Digital 

223 Nanodot Gradients (DNGs) using a unit cell approach. (A) Schematic of a linear 

224 DNG and unit cell parameters for the i th column of nanodots. (B) Calculation of unit cell 

225 dimensions at each position through an iterative process using the cumulative sum of 

226 unit cell dimensions to determine position along the length. 
227 
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228 The width of the gradient is divided by d, to estimate the number of nanodots in the 

229 column. The value of d, is then recalculated to equally space the integer number of 

230 nanodots along the width, ensuring that the distribution of nanodots is symmetrical. The 

231 unit cells are concatenated to form a column of points in the y direction with constant 

232 spacing as multiples of of,-. The x coordinate for each column is the position / at which the 

233 unit cell dimensions are calculated. Thus, using this algorithm it is easy to form a 

234 gradient with any slope as the size of the unit cell is directly derived from the value of the 

235 density function at a specific position. Fig. 2 shows a linear density curve spanning from 

236 a minimum density D min =0.01 to a maximum density D max =0.30 (Eq. 2) with 

237 corresponding unit cell dimensions along the length of the gradient. Eq. 2 is normalized 

238 to span the dynamic range over the total length L=100 /jm of the gradient. 
239 

240 D(l) = (D -D . )— + D . (2) 

\" / V max min / mm \ I 

241 
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Figure 2: Programmed and measured density for an ordered linear DNG. A 100 * 

100 jjm 2 linear DNG comprised of 0.04 pm 2 nanodots and spanning from 0.01 to 0.30 
density. The measured density (left axis) precisely matches the programmed density 
(R 2 =0.9998). The unit cell dimensions (right axis) change with the inverse square root of 
the density. Top insets show various spacing of nanodots at positions along the DNG. 



To verify that ordered gradients matched the programmed density functions, gradients 
were exported as bitmap images and the ratio of black to white pixels, averaged over 
several columns of nanodots, was used to measure density. The resulting gradients 
were found to match the programmed density functions with high fidelity (R 2 =0.9998). 
Using this algorithm, a density of zero cannot be reached, as it would require an infinitely 
large unit cell. Consequently, when designing gradients of fixed length, the lowest 
density that is accessible may be limited since large unit cells are required at very low 
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256 density. 
257 

258 Random Gradient Algorithm 

259 Naturally occurring gradients appear continuous at the microscale, but are in fact digital 

260 at the scale of molecules and proteins. The diffusion of biomolecules is subject to 

261 Brownian motion, and is therefore expected to be random. In an effort to mimic the 

262 nanoscale noise present in vivo, we present an approach to produce DNGs with fully 

263 randomized nanodot positions. Similar to ordered DNGs, random gradients have 

264 increasing density along the length and constant density along the width. 
265 

266 Several strategies were considered to produce randomized DNGs. The first was akin to 

267 the one proposed by Dalby, ef al. [1 6] and consisted of starting with a regular distribution 

268 of dots, and randomly displacing the nanodots within their "unit cell" so as to avoid 

269 overlap with neighboring dots. However because this approach only produces a quasi- 

270 random distribution, it was not pursued. A second approach that was considered was to 

271 start with a pseudo-random dot pattern, and then compensate for the overlap by locally 

272 moving dots. This parallels a strategy used in back-lit displays [29], however these 

273 patterns are unsuitable for our application due to the dependence on or optimization for 

274 quasi-random distributions. Whereas we don't see fundamental obstacles in developing 

275 an algorithm to randomly redistribute overlapping dots, such algorithms tend to be 

276 computationally intensive as the position of each dot needs to be tested against each 

277 other dot, and require multiple iterations to achieve the final distribution [27]. We thus 

278 devised a novel algorithm to create pseudo-random nanodot patterns for cell navigation 

279 that could translate arbitrary density functions into random nanodot patterns, and that 

280 was computationally economical. 
281 
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282 The algorithm we developed compensates for the random overlap of dots by adding an 

283 excess of dots to achieve the programmed density. The algorithm first defines equally 

284 sized boxes to smoothly replicate the programmed density function without visible steps. 

285 The number of nanodots within each box was calculated based on the density function 

286 with predicted overlap at the box's position, and random x and y coordinates were 

287 generated for each nanodot. Only the nanodot coordinates were stored. Noteworthy, the 

288 shape of the nanodots can be defined subsequently as the exported CIF file allows 

289 drawing shapes such as circles, polygons, etc. at the stored coordinates, thus making 

290 this approach flexible while reducing the file size and computation time. A consequence 

291 of the random nanodot position is that overlap between nanodots is possible. The 

292 frequency of overlap increases with density, which in turn results in a lower density than 

293 expected based on the surface area of all nanodots. However, because the distribution 

294 of the nanodots is random and the density known, the overlap is predictable, and can be 

295 compensated for, so that the designed DNG matches the density function with high 

296 fidelity. For each box with area A box , we first compute the probability P cov that any given 

297 point in the box is covered by a nanodot of area A dot (Eq. 3). 
298 

299 P C ov = y^ ( 3 ) 

"■box 

300 

301 The probability P not cov that this point will not be covered by the nanodot is simply 

302 1-Poov (Eq.4). 
303 



304 P 



305 



f A A A -A 

"■box "dot 



V "box i 



"box 



(4) 
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306 The probability that this point will not be covered by N nanodots simultaneously can then 

307 be found (Eq. 5). 
308 



309 P t 

not cov 



320 



f A -A A 

"-box "-dot 
V \ox J 



(5) 



310 

311 To determine the total area of the box covered by nanodots (A cov ), the probability that a 

312 point will be covered is integrated over the area of the box (Eq. 6). 
313 

314 A =[ A ""'(1-P )dA (6) 

cov Jq V not cov/ > ' 

315 

316 A cm = A h0X - (Ab '"~ A f jN (7) 

A box 

317 

318 Eq. 7 can then be solved for N to determine the number of nanodots required for a given 

319 A cov 'm each box. 



log(l-A cov /A fcojc ) 



321 N — — — — — (8) 

\og{l-A d JA box ) 

322 

323 Using Eq. 8 it is possible to obtain the number of nanodots required for any given density 

324 by substituting D for A cov /A box . As the area occupied approaches the area of the box, the 

325 number of nanodots required increases rapidly, and is infinite for a density of 1. The 

326 maximum density for a DNG was thus set to 0.9999. Without compensation, at density of 

327 0.9999 would require 9,999 nanodots of 200 * 200 nm 2 to be seeded into a 400 x 1 urn 2 
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328 box to reach this density. When accounting for the overlap, 92,099 nanodots are 

329 required, roughly a tenfold increase. Thus, producing a typical linear gradient, ranging in 

330 density from 0.0002 to 0.4444 requires 1,061,307 nanodots and a computation time of 

331 4.56 min. The ordered algorithm for the same linear gradient requires 885,737 nanodots 

332 and a computation time of 3.88 min. Therefore, the random gradient algorithm overall 

333 required a 19.8% increase in nanodots and a 17.5% increase in computation time using 

334 our workstation. It is expected that the time will increase for random DNGs that reach a 

335 higher density. 
336 

337 We compared the programmed density with the measured density of a monotonic 

338 exponential curve with decay constant k (Eq. 9) using both the ordered and random 

339 gradient generation approaches. Eq. 9 is normalized to span the dynamic range over the 

340 length of the gradient. The random gradients produced with this algorithm were found to 

341 follow the programmed curve with high fidelity (R 2 = 0.99986), and accurately match the 

342 measured density of ordered gradients produced with the same input function (Fig. 3). 
343 

e L -1 

344 D(l) = (D -D ) — + D (9) 

VV V^max mm/ k 1 mm V ' 

e - 1 
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Figure 3: Ordered and random DNGs superposed with their programmed 
exponential density function and actual, measured density. A 100 x 100 pm 2 

exponential (k=3) DNG spanning from a density of 0.01 to 0.30 is shown in ordered (top 
half) and random (bottom half) form. The random DNG was programmed by subdividing 
it into 1 urn wide boxes with random seeding of nanodots and compensating for overlap. 
The measured density follows the programmed exponential curve with high fidelity for 
both ordered (R 2 =0.99996) and random (R 2 =0.99986) DNGs. 



The distribution of neighboring nanodots can be parsed to verify whether it satisfies the 
conditions of randomness. Since density changes along the length of the DNG, 
randomness can only be assessed in the direction of constant density, perpendicular to 
the gradient. The randomness along each box was verified using Ripley's K function 
based on the number of nanodots N pi within a distance s from the /* point p, taken over 
the sum of all points n and normalized by the area A (Eq. 1 0). 



361 



K(s) = -J j N p is)/X 



n z 



(10) 
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362 

363 For a homogenous, random Poisson distribution, the expected value of Eq. 10 is tts 2 . 

364 Deviations from tts 2 indicate regions of clustering or dispersion [31]. Using the 

365 coordinates of nanodots from the DNGs, K(s) was found to lie within a 95% confidence 

366 interval obtained from 10 simulations of randomly distributed coordinates at the same 

367 density. This suggests the nanodots are spatially homogenous and randomly distributed 

368 compared to ordered gradients, which lie outside the confidence interval (Fig. 4). 
369 

20 
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5 
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370 

371 Figure 4: Ripley's K analysis for randomness and spatial homogeneity of a sample 

372 box produced with the ordered and random DNG algorithms. K(s) for ordered and 

373 random gradients for a 10 * 400 pm 2 box at 0.20 density are shown. A 95% confidence 

374 interval was calculated as 1.96 times the standard deviation of K(s) from 10 simulations 

375 of randomly distributed points at 0.20 density. The inset shows a magnified portion of the 

376 graph indicated by the box. 
377 

378 Non-monotonic gradients 

379 We have shown that linear and exponential gradients can be produced using either the 

380 ordered or random algorithms. These curves are monotonic, meaning they only ever 

381 increase or decrease. Given the approaches outlined here, more complex gradients can 
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382 be easily generated from any input density function, specifically non-monotonous curves. 

383 To study the ability of cells to recognize an average gradient in a non-monotonous 

384 environment, we propose linear and exponential gradients that are superposed with 

385 sinusoidal functions with either constant amplitude, or with linear or exponentially 

386 changing amplitude. In the anticipation of future cell experiments, we also programmed a 

387 series of sinusoids with varying frequency and amplitude, while having an average slope 

388 of zero to act as negative controls. 
389 

390 Here, we demonstrate the flexibility of the algorithms with the production of non- 
391 monotonic gradients. One such gradient produced is a sinusoid with exponentially 

392 increasing amplitude superposed an exponential gradient. Eq. 11 gives the input density 

393 function for such a curve, where A is the amplitude, 6 is the number of oscillations, and 

394 k-i and k 2 are the decay constants for the average gradient and amplitude of the sinusoid, 

395 respectively. 
396 
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399 

400 Figure 5. Non-monotonic random DNG superposed with its density function and 

401 actual, measured density. A 100 x 100 pm 2 non-monotonic DNG. The programmed 

402 density function (line) is a sinusoid (A=0. 10, B=20) with an average exponential trend 

403 (ki=3) and exponentially increasing amplitude (k 2 =3). Measured density from bitmap 

404 (circles) accurately follows the programmed function (R 2 =0.9988). Insets show close-up 

405 views of nanodots and reveal increasing overlap of nanodots at higher densities. 
406 

407 One-Hundred-Gradient Array 

408 The flexibility of the gradient algorithms and the fabrication method discussed below was 

409 leveraged by producing an array of 100 distinct gradients within a 35 mm 2 area (Fig. 6, 

410 Table. S1). The array is comprised of 20 ordered and 80 random gradients, with 

411 densities ranging from 0.0002 to 0.4444 for a maximum dynamic range of 3.85 OM. For 

412 ordered gradients, this corresponds to a maximum pitch of 14.8 pm at low density, and 
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413 100 nm at high density. The minimum density is limited by the size of the cell and its 

414 capacity to sense the gradient, i.e. the cell will not sense the gradient if dot spacing 

415 exceeds one third the size of the cell, which we previously assessed using the migratory 

416 response of C2C12 myoblasts [14]. The maximum density is defined by the resolution of 

417 e-beam lithography, which limits the minimum pitch to 100 nm. To address how 

418 randomness affects cell navigation, 20 gradients (10 linear and 10 exponential) were 

419 produced as both ordered (#1-20) and random (#81-100) gradients. This portion of the 

420 array will serve to address gradient sensing mechanics for cell migration that may arise 

421 from either (i) the absolute concentration of the gradient at a given location or (ii) the 

422 concentration ratio between the cell's leading and trailing edges. 
423 

424 The other 60 gradients are non-monotonic and random. These consist of 14 sinusoids 

425 with no average slope (#21-34) that serve as controls. These vary in frequency and 

426 amplitude, and may have constant or increasing amplitude along their length. Controls 

427 for ordered gradients were designed in our prior work. The remaining 46 gradients (#35- 

428 80) have sinusoidal curves with various levels of complexity. These include sinusoids 

429 with linearly (#41-60) and exponentially (#61-80) increasing average density, and non- 
430 monotonic functions with linearly increasing amplitude and average density (#35-36) and 

431 exponentially increasing amplitude and average density (#37-40) to demonstrate the 

432 flexibility of this new approach in gradient generation. As discussed, the sinusoidal 

433 curves with different amplitude and frequency may introduce obstacles for cells trying to 

434 sense the overarching gradient, and replicate effects that may occur due to cell 

435 "mosaicism" in a controlled manner. 
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1 Gradient Array Fabrication 

2 The hundred-gradient array was etched 100 nm deep into a Si wafer by e-beam 

3 lithography (Fig. 7). The integrity of individual dots for ordered and random gradients 

4 was confirmed by scanning electron microscopy (SEM) (Fig. 8). 
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Figure 7: Optical images of the 1 00-gradient-array. The 100 DNGs were patterned 
into a Si wafer using e-beam lithography. (A) Image of all gradients; scale bar is 1 mm. 
(B) Dark-field image of DNG 67 which is a random sinusoid with exponentially increasing 
average density; scale bar is 100 pm. 
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12 Figure 8: Nanodot distribution in ordered and random DNGs. SEMs of ordered (A-D) 

13 and random (E-H) exponential gradients at low (A,B,E,F) and high (C,D,G,H) densities. 

14 (H) highlights the random overlap between nanodots that is compensated for by the 

15 random DNG algorithm. Scale bars are 10 urn (top row) and 500 nm (bottom row). 
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16 To translate the etched Si wafer into substrate-bound protein gradients, we employed 

17 lift-off nanocontact printing [14]. First, a PDMS intermediate replica was produced, 

18 followed by a second replication into low-cost Norland Optical Adhesive. A flat PDMS 

19 stamp was inked with a protein solution, and through contact with the plasma-activated 

20 NOA lift-off stamp, a monolayer of protein was selectively removed from the surface of 

21 the flat stamp leaving behind the digital nanodot protein pattern. Next, the flat PDMS 

22 stamp was brought into contact with a plasma-activated glass slide to transfer the DNG 

23 pattern. To confirm the accuracy of the replication and printing process, images of the 

24 design and the printed protein DNG were digitally overlapped and compared (Fig. 9), 

25 revealing a high fidelity between the two patterns. 
26 
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28 Figure 9: The printed nanodot pattern accurately replicates the design. Bitmap 

29 image of the programmed design colored in red (A, D) compared with a fluorescent 

30 image of nanocontact printed IgG by lift-off (B, E, green) and merged (C, F). In the inset 

31 (E), out of -1000 spots, 21 are missing, indicating that the replication works well. Scale 

32 bars A-C are 10 pm, scale bars D-F are 2 pm. 
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33 

34 The pattern overlay in Fig. 9 was characterized by thresholding the printed image with 

35 boundaries of 31 and 255 in ImageJ and comparing it to the bitmap image (Fig. S3). The 

36 total number of nanodots is difficult to determine owing to the random overlap, but we 

37 estimate that there are -1000 spots, and using the above threshold, 21 were not found 

38 in the print. Overall, the transfer process is thus accurate to -98% in this example. Dust 

39 particles on the Si wafer or on the intermediate PDMS replica, or mechanical damage 

40 due to the replication process could account for the absence of these protein dots. 

41 Whereas the fidelity of the replication and printing process is thus high, it might be 

42 improved further by employing more durable polymers during the replication process as 

43 well as by working in a cleanroom environment throughout. 
44 

45 Conclusion 

46 Patterned substrate-bound protein gradients are a valuable tool to study a number of 

47 biological processes such as neuronal development or regeneration. The novel 

48 algorithms presented here, providing pseudo-randomness at a scale commensurate to 

49 the cellular level, can provide control over geometry and noise in DNGs. Combined with 

50 high-throughput patterning, an array of one hundred DNGs with linear, exponential, and 

51 non-monotonic slopes featuring 57 million spots over an area of 35 mm 2 can be 

52 patterned at once in a matter of minutes. The diversity of DNGs shown here will help 

53 study and quantify the mechanisms by which cells sense and navigate through 

54 immobilized gradients. There are many opportunities for refining digital nanodot patterns. 

55 Firstly, to mimic the self-repellent nature of proteins adsorbing to surfaces, it might be 

56 useful to develop an algorithm that prevents, or limits, the overlap of nanodots. 

57 Secondly, whereas here two sinusoidal curves were superposed in one direction, such 

58 curves might be generated in two perpendicular directions and create two- dimensional 
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59 navigation landscape while better mimicking local accumulation of guidance cues as 

60 puncta. Thirdly, it might be possible to program density functions that introduce clusters 

61 of noise that more accurately replicate the noise and mosaicism imposed by individual 

62 cells in vivo which can extend over tens of micrometers. Fourthly, it should be possible to 

63 pattern overlapping DNGs of different proteins that run in the same, or different 

64 directions, as well as produce any type of navigational landscapes that are found in vivo 

65 simply by converting the recorded densities into digital nanodot patterns, following the 

66 trend of rapid prototyping of replicas of living tissues [32]. Before expanding the nanodot 

67 patterning, it will be important to test and validate the current DNGs and establish the 

68 optimal conditions along with the suitable non-patterned reference surfaces for each 

69 study [33]. 
70 
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ABOVE THRESHOLD 
(SET TOO- WHITE) 



BELOW THRESHOLD 
(SET TO 1- BLACK) 



Figure S1: Random gradients produced by a random matrix threshold approach. A 

matrix of pseudo-random numbers is generated. Values greater than or equal to the 
density threshold are set to 0 (white), while values less than the threshold are set to 1 
(black). The binary array can then be exported directly as a bitmap image file. Each 
nanodot is represented by one pixel. Thus, for a 400 * 400 pm 2 sized area with 200 x 
200 nm 2 nanodots, a 2000 x 2000 matrix with 4 million values is required. This approach 
does not provide a fully random configuration since nanodots are aligned to a grid. While 
these patterns appear random to the eye, the underlying grid might be sensed at the 
cellular scale. Theoretically, it is possible to further randomize the position of dots by 
subdividing the area, e.g. using a 50 nm grid to position and draw 200 x 200 nm 2 pixels, 
but this would come at the cost of increased computational requirements and the 
possibility of overlap between adjacent dots. 
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Figure S2: Ordered and random measured density matches programmed density 
overfull range. Linear gradients from densities of 0.01 to 0.99 were shown to match 
(dots) the programed functions (line) with high fidelity for both ordered (R 2 = 0.9988) and 
random (R 2 = 1.0000) gradients. Using either approach, a high dynamic range can be 
achieved with near perfect match to the programmed function. 
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Figure S3: Image processing procedure to assess alignment of the DNG design 
and print. The fluorescent image of nanocontact printed IgG was first thresholded in 
I mage J with boundary values of 31 and 255. The image was then transformed to binary 
and the binary values inverted to facilitate visualization. The edited fluorescent image 
(green) was then merged with the bitmap (red) and yellow dots, indicative of non-printed 
dots, were counted to determine to what extent the print matched the design. 
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